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ABSTRACT 

We have searched the WISE first data release for widely separated (^10,000AU) late 
T dwarf companions to Hipparcos and Gliese stars. We have discovered a new binary 
system containing a K-band suppressed T8p dwarf WISEP J1423+0116 and the mildly 
metal poor ([Fe/H]=-0.38±0.06) primary BD+01 2920 (Hip 70319), a Gl dwarf at a 
distance of 17.2pc. This new benchmark has T o g=680±55K and a mass of 20— 50Mj up . 
Its spectral properties are well modelled except for known discrepancies in the Y and 
K bands. Based on the well determined metallicity of its companion, the properties of 
BD+01 2920B imply that the currently known T dwarfs are dominated by young low- 
mass objects. We also present an accurate proper motion for the T8.5 dwarf WISEP 
J075003.84+272544.8. 

Key words: surveys - stars: low-mass, brown dwarfs 



1 INTRODUCTION 

An accurate understanding of the physics of ultra-cool atmo- 
spheres (T e ff < 2300K) is a m ajor and ongoing challenge for 
theo ry (e.g . Allard ct al. 1997). Complex molecular opacities 
(e.g. iBarber et ai1l2006lh cond ensate clouds and their prop- 
erties (e.g. Allard et al.| [200D. and n on-equilibrium chem - 
istry (i.e. vertical transport or mixing; ISaumon et al.|[2007l ) 
are significant sources of uncertainty in the models. How- 
ever, it is crucial to improve our understanding if we are to 
effectively measure the properties of substellar populations 
(brown dwarfs and giant planets) and study their formation 
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and evolution (e.g. iBate et al.l 


2002:lGoodwin & Whitworthl 


2007: Stamatcllos et al. 2007: 


Sumi et alj 201lh. 



Building on the samp les of L (~2300K -1500K) and 
T dwarfs (~1500K-500K) (|Kirkpatrickl 120051) identifi ed in 
the Two- Micron All Sky Survey dSkrutskie et al-I^OPd ), the 
DENIS survey dEpchtein et all Il997f ) and the Sloan Dig- 
ital Sky Survey ( York et al.l l2000t ). a new generation of 
infrared surveys is expanding our search-space into much 
greater vo lumes. The UKIRT I nfrared Deep Sky Survey 
(UKIDSS; lLawrence et "all 120071 ) is sensitive to mid L - 
mid T dwarfs out to ~100pc over ~15% of the sky. The 
VISTA surveys (e.g. the Viking and VHS surveys) will ex- 
pand this coverage to ~50% of the s ky in the next few 
years. For the latest T dwarfs (T8-9; e.g.lWarren et al.ll2007l: 
iBurningham et ail 120081 : iDelorme et al. 20081 : Lucas et al.l 
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120101 ; iBurningham et all l2011ch the sensitivities of these 
surveys are matched by those of the Wide- field Infrared 
Survey Explorer (WISE; IWright et all l201p| ). probing to 
distances of ~15-25pc. And for even lower temperatures 
(T c s <500K) an increased mid-/near-infrared flux ratio 
is allowing WISE to uncover the new Y dwarf class in 
the ~300-500K range (|Cushing et al.ll201lh . Together these 
surveys are characterising a rapidly growing population 
in the near- and mid-infrared (e.g. lLodieu et al.l 120071: 
Pinfield et all 120081; IBurningham et al.l l2010bl ; iRevle et al] 



20ld ; iKirkpatrick et alj|201lj ) 



With sensitivity to a growing field L-, T- and Y- 
dwarf search-space i t is becoming feasible to s e arch for mul- 
tiple s ystems (e.g. [B urningh am et al. 2009a ; Zhang et al] 
20ld: 'IBurningham et al.l koiOal: iDav- Jones et al.l l2oilT 



Murray et al.ll201ll : IBurningham et al.ll2011al: iLeggett et all 
2010d') or moving group a ssociations fe.g. lClarke et al.ll20ld ; 
Galvez-Ortiz et all l20ld ). The physical properties (mass, 



age and metallicity) of such objects can be constrained 
through association with more readily characterisable stel- 
lar companions or moving group members, establishing them 
as benchmark objects that can test the theory or more di- 
rectl y map physical pro perties onto spectral characteristics 
fe.g. lPinfield et al.ll2 006l. Indeed, by searching for even rarer 
benchmarks with better physical constraints which span a 
more extreme range of properties, it will be possible to pro- 
vide the strongest tests for the model atmospheres, a goal 
that absolutely requires sensitivity to large volumes. 

In this paper we present a search of the WISE first data 
release for widely separated late T dwarf companions to stars 
with known parallaxes. Section 2 describes our WISE sample 
selection, and Section 3 the method used to identify candi- 
date binary associations. Sections 4 and 5 present our spec- 
troscopic and additional photometric data, and in Section 6 
we derive candidate proper motions. Section 7 statistically 
assesses the expected level of false positives in our search, 
and Sections 8, 9 and 10 discuss the properties and charac- 
teristics of a newly discovered benchmark system. Conclu- 
sions and future work are in Section 11. 



2 WISE SAMPLE 

We identified candidate mid-T and later type objects in the 
WISE Preliminary Data Release source catalogue, which we 
accessed via the NASA/IPAC Infrared Science Archive's 
catalogue query engine. We performed a series of all-sky 
searches using structured query language input to select 
sources with constraints on signal-to-noise and colour, and 
with detections in various combinations of bands chosen to 
optimise sensitivity to late T dwarfs. We always required 
a detection in the W2-band with signal-to-noise (SNR) of 
at least 10. If Wl- and W2-band detections are available 
then we require Wl — W2 ^2.0 to select spectral ty pes 
of >T5 l|Kirkpatrick et aT. 1 12011| ; iMainzer et al] l201lf ). If 
W2- and W3-band detections are available then we require 
W2 — W3 ^2.5 in order to avoid du sty galaxies such a s 
ULIRGS, LINERS and obscured AGN (|Wright et al.ll20ld ). 
As well as our WISE-band detection requirements we also 
made use of the WISE catalogue cross-match with the 
2MASS point source catalogue, to divide our searches into 
objects that are detected in 2MASS (within 3 arcseconds 



of the WISE position) and those that are not. For 2MASS 
detected objects we required that either H — W2 J?2.5 or 
J — W2 ^3.5 so as to remove L and early T dwarfs. The full 
complement of searches and the number of sources identified 
in each is shown in Table [1] 

The search requiring detection in only the W^-band will 
be the most sensitive to faint objects with red WISE colours 
since the WISE sensitivity limits (all-sky 5- u Vega limits 
are 1 ^1=16.5, W2=15.5, W3=11.2, W4=7.9; IWright et all 
mean that objects with Wl — W2 >2 will generally 
only be detected in W2 for W2=14.5-15.5 (i.e. at least ~75% 
of the W2 survey volume). The other multi-band combina- 
tions cover the full range of detection/non-detection combi- 
nations that might be expected for T dwarfs. 

For c omparison, the re c ent la rge-scale WISE search 
made by IKirkpatrick et all (120111 ) overlaps significantly 
with our search-space. However, they use a slightly bluer 
Wl — W2 selection, and where we require SNR>10 

in the W2-band they require at least 8 separate detections 
(SNR>3) in the individual W2 exposures. 



3 IDENTIFYING CANDIDATE BINARY 
SYSTEMS 

To identify candidate binary systems we cross-matched our 
candidate late-T sample with a list of potential primary 
stars with measured parallax distances, and imposed sepa- 
ration constraints on the potential binary pairings as well as 
absolute magnitude constraints on the candidate T dwarfs 
(where we assumed a common distance for components). 
The list of potential primary stars was mad e by combin- 
ing t ogether the latest Hipparcos catalogue l|van Leeuwenl 
l2007h with the most recent ve rsion of the catalogue of nearby 
stars l|Gliese fc Jahreisslll99ll) . Hipparcos provides astromet- 
ric measurements (in position, parallax, and annual proper 
motion) with uncertainties in the range 0.7-0.9 milliarcsec 
(mas) for stars brighter than V=9. The catalogue is on the 
ICRS reference system and has proper motions consistent 
with an inertial system at the level of ±0.25mas/yr. The 
Third Catalogue of Nearby Stars (CNS3) contains informa- 
tion on all known stars within 25 parsecs based on an ex- 
tensive literature search during almost four decades. 

We required on-sky angular separation of candidate 
pairs to be ^300 arcseconds to reduce contamination from 
random alignments, and physical separation of ^10,000AU 
(where we use the distance of the primary to convert angular- 
separation into a physical line-of-sight separation) , since the 
great majority of known wide ultracool stellar companions 
have separations b elow this limit (e.g. IZhang et al.l l20ld ; 
iFahertv et"ail l2010). In addition we use the known distance 
of each candidate primary to estimate the Mwi that the 
T dwarf candidate would have at this distance, and reject 
any associations where the T dwarf candidate would have 
Mw 2 ^11.5 (i.e. targeting late T dwarf companions, e.g. Fig 

4 of ILeggett et al] [2010a . where M[ 4 5 ] is a good proxy for 
Mwi)- The candidate T dwarf components of the possible 
binary systems are distributed within our WISE source se- 
lections as summarised in Table[T] These T dwarf candidates 
were visually inspected using the WISE image server at the 
NASA/IPAC Infrared Science Archive, and candidates re- 
jected if the source did not appear point-like in any of the 
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a W2 signal-to-noise always >10 (w2snr>=10). 

b Non-detections denned as (w*mpro is null or w*sigmpro is null). 

c Detections defined as (w*mpro is not null and w*sigmpro is not null). 

d 2MASS non detections (tmass_key is null). 

e The numbers in brackets are for candidates that passed visual inspection. 
1 2MASS detections with H — W2 >2.5 or J - W2 >3.5. 

Table 1. WISE late T candidate sample. Twelve separate selections were made, six requiring non-detection in 2MASS and six requiring 
2MASS detection with red 2MASS-WISE colour. The number of candidate T dwarfs, and those that became wide companion candidates 
(see Section 3) are indicated for each search, where various combinations of detection and non-detection were explored in the four WISE 
bands. 



bands, formed part of a blended structure, or was clearly an 
artefact (e.g. part of a diffraction spike). 

Five candidate binary systems passed visual inspection. 
One W2-only detected candidate remains an unconfirmed 
interesting candidate without any additional survey data 
(e.g. UKIDSS, VISTA) to facilitate proper motion measure- 
ments. The other 4 are listed below: 

• WISEP J075003. 84+272544. 8 is a W2-only detected 
candidate 265 arcseconds from the star Hip 38228, a G5IV 
star at 22pc. This can didate is a known T8.5 d warf discov- 
ered (with WISE) by iKirkpatrick et al.l (|201lh . though we 
subsequently show (Section [©J that it is not a companion to 
Hip 38228. 

• WISEP J142320.86+011638.1 (WISEP J1423+0116) is 
a W^2-only detected candidate 153 arcseconds from the star 
Hip 70319 ( BD+01 2920), a G 1 V sta r at 17.2pc. It was not 
identified by IKirkpatrick et al.l (|20 1 it) because it is only de- 
tected in 7 separate W2 exposures in the WISE Preliminary 
Data Release. This T dwarf is the main subject of this paper. 

• WISEP J145715.85-212207.6 is a Wl + W2 + W3 de- 
tected candidate near the system Gl 570ABC (Hip 73182 
and Hip 73184), a K4V+M1.5V+M3V triple system. This 
candidate is a known ( discovered in 2MASS) T8 member 
of the multiple system (|Burgasser et al.ll200oh . The WISE 
catalogue does not list the source as a 2MASS detection be- 
cause its high proper motion leads to the WISE and 2MASS 
positions being separated by more than 3 arcseconds. 

• WISEP J150457.58+053800.1 is a Wl + W2 detected 
candidate 63 arcseconds from Hip 73786 (GJ 576), a K8V 
star at 18.6pc. This candidate is a known (discovered in 
UKIDSS) T6p companion to this somewhat metal poor star 
i|Murrav et al.ll201ll ; IScholdbOlOr ,. 

Figure [T] shows J- (VISTA) and VF2-band images for 
WISEP J1423+0116, and indicates its separation from the 
nearby high proper motion star Hip 70319 (BD+01 2920). 



4 SPECTROSCOPY 

Near-infrared spectroscopy of WISEP J1423+0116 
(BD+012920B; see Section [7J) was obtained using the 
Gem ini Near InfraRed Spectrograph (GNIRS; lElias et al.l 
2006) mounted on the Gemini-North telescope on the 
night of 16 th May 2011. The target was observed in 
cross-dispersed mode capturing the full 0.8-2.5u.m region 
with a 1.0" slit delivering a resolving power of R~ 500. The 
data were redu ced using GNIRS routi nes in the Gemini 
IRAF package (|Cooke fe Rodgersll2005h . using the nearby 
F5V star Hip 63976 for telluric correction. The telluric 
standard spectrum was divided by a black-body spectrum 
of an appropriate T e a after removing hydrogen lines by 
interpolating the local continuum. The rectified standard 
spectrum was then used to correct for telluric absorption 
and to provide relative flux calibration. The overlap regions 
between the orders in the V— ,J— and //—bands agreed 
well suggesting that the relative flux of the orders is well 
calibrated. The resulting YJHK spectra are shown in 
Figured 

In Figure [2] we compare our GNIRS spectrum of the new 
T8 with that of the T8 spectral templa te 2MASS J04151954- 
0935066 from iBurgasser et alj (|2006bl ). The close similarity 
of the spectra over most of the wavelength range argues 
strongly for T8 classification, which is reflected in the val- 
ues found for the spectral typing flux ratios (see Table [2}. 
Although the new T8 closely traces the T8 template over the 
0.9 — 1.9(J.m range, it displays a considerably more depressed 
AT-band flux, which is interpreted as due t o strong colli- 
siona lly induced absorption by H2 (CIA H2; ISaumon et alj 
Il994h . Increased CIA H2 is typically attributed to higher 
pressure atmosph eres arising from lower-metallicity and /or 
high-gravity (e .g. IBurgasser et al.|[2002l ; iKnapp et al. I l2004l ; 
iLiu et ali r2007). For this reason we assign the type T8p to 
WISEP J1423+0116, where the 'p' suffix denotes it is pecu- 
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WISE Channel 2 

Figure 1. J-band and W2-band images of WISEP J1423+0116. 
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Figure 2. GNIRS spectra for WI SEP J1423+0116 (BD+ 012920B; see Section compared to T8 spectral type template 
2MASS J04151954-0935066 taken from lBurgasser et alj j2006rj) . The error spectrum is shown offset by -0.1. 



liar, alluding to the poor match with the template in the K 
band. 



5 NEW PHOTOMETRY 

Database photometry of WISEP J1423+0116 was obtained 
from the WISE Preliminary Data Release catalogue, the 
WFCAM Science Archive (UKIDSS Large Area Survey) and 
the VISTA Science Archive (VIKING proprietary data ac- 
cess). In addition, observations were taken at the Telescopio 
Nazionale Galileo (tng) and with the Spitzer Space Tele- 
scope in its warm phase. 



Near- infrared photometry was me asured using the N ear 
Infrared Camera Spectrometer (NICS; iBaffa et alj|200 if ) at 
the 3.58-m optical/infrared TNG located on La Palma, on 
the night of the 7 th of May 2011 for the H band and the 
night of the 10 th of June 2011 for the Y band. The data 
were obtained in large field mode, with a pixel scale of 0.25 
arcsec/pixels and a field of view of 4.2x4.2 arcmins. The 
data were processed using the NICS science pipeline SNAP 
provided by the TNG. /f-band observations consisted of a 
50 point jitter pattern with individual 10s exposures and 6 
co-adds per jitter point, accumulating to a total exposure 
time of 50 minutes. In the Y-band a 10 point jitter was used 
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Table 2. The spectral flux ratios for WISEP J1423±0116 
(BD+012920B; see Section [TJl. The locations of the numerators 
and denominators are indicated on Figure [2] 



for the same exposure time and co-adds, resulting in a total 
exposure time of 10 minutes. We calibrated each image onto 
the MKO system using ~30 field stars in the frame. 

Warm-Spitzer photometric data were obtained for 
WISEP J1423+0116 on the 21 st August 2011, via the Cy- 
cle 7 GO program 70058. Individual frame times were 30 s 
repeated six times, with a 16-position spiral dither pattern, 
for a total integration time of 48 min in each of the [3.6] and 
[4.5] bands. The post-basic-calibrated-data mosaics gener- 
ated by version 18.18.0 of the Spitzer pipeline were used to 
obtain aperture photometry. The photometry was derived 
using a 7-arcsec aperture and the aperture correction was 
taken from the IRAC handbook. The error is estimated by 
the larger of either the variation with the sky aperture or 
the error implied by the uncertainty images. 

Tables [3] and U contain the available photometry and 
colours, respectively, for the T dwarf. We present WISE pho- 
tometry where the signal-to-noise is positive and note that 
the Wl and W3 magnitudes are brightness upper limits. The 
near infr ared photometry is on the Mauna-Kea Observatory 
system l|Leggett et al.l 120061 ) except for the TNG y-filter, 
which is slightly different (A c =1.02urn, FWHM=0.13um) to 
the MKO Y filter (A c =1.02um, FWHM=0.10um). In the 
absence of a measured if -band magnitude we have deter- 
mined a synthetic J — K colour using our GNIR S spectra 
and a spectrum of Vega (Bohlin & Gilliland 2004) both con- 
volved with the re sponse functions for the pass-bands (e.g. 
iHewett et al . 2006). This synthetic colour (see Table© com- 
bined with the J-band magnitude produced our K band es- 
timate. For the mid infrared photometry we note that while 
similar to Wl and W2, the Spitzer [3.6] a nd [4.5] bands 
have some significant differences (see Fig 2 of iMainzer et al.l 



l201ll ). There are multiple measurements of Y-, J- and H- 
band photometry, though no evidence of variability is seen 
(to within the uncertainties) in the photometric brightness. 



6 PROPER MOTIONS 

We measured the proper motion of WISEP J1423+0116 us- 
ing a VISTA VIKING image of April 2010 and two UKIDSS 
images from May 2008 (with lower signal-to-noise of ~6). 
This avoids using the larger point-spread-functions inher- 
ent in the WISE images (~6.5 arcsecs in W2). The base- 
line between the two near infrared epochs was 1.89 years. 
We took the measured x,y coordinates from the standard 
CASU pipeline reductions of all images and using 59 objects 
within 4 arcminutes of the target, transformed the UKIDSS 
frames onto the standard coordinate system of the VIKING 
frame using a simple linear model. The relative proper mo- 
tion for all objects were found from linear fits to the stan- 
dard coordinates at the different epochs. A correction to 
an absolute system was estimated from the median differ- 
ence between measured relative proper motions and 6 SDSS 
objects in the fiel d with proper motions in the catalog of 
iMunn et all (|2004h ■ The derived proper motion for WISEP 
J1423±0116 was corrected for an assumed parallax of 50mas 
(see Section [7J), and final uncertainties are based on the 
formal uncertainties of the measured coordinates combined 
with an additional allowance for the centroiding accuracy in 
the low signal-to-noise LAS image (~ ±0.5 pixels estimated 
using Monte Carlo techniques) leading to a proper motion 
uncertainty of ±50 mas/yr. The proper motion of WISEP 
J1423±0116 is fi acosS = 261±56 mas yr" 1 , fig = -444±52 
mas yr -1 , which is within 0.7a of the Hipparcos proper mo- 
tion vector of Hip 70319 (BD+01 2920; ^ QC os«5 = 223.8 + 0.4 
mas yr -1 , fi 6 — —477.4 ± 0.4 mas yr -1 ). These objects are 
thus a common proper motion pair. 

We also measured the proper motion of WISEP 
J075003. 84+272544.8 using two UKIDSS LAS J-band 
epochs with a baseline of 2 years. We applied a second or- 
der polynomial transformation between the two epoch im- 
ages to correct for any non-uniformity in the focal plane. 
Seventeen reference stars (J <18.1) were used, distributed 
around the target (with at least 3 per quadrant) with sepa- 
rations within 2 arcminutes. A correction was applied to an 
absolute system using apparent proper motions of nearby 
galaxies. The uncertainties were calculated using the stan- 
dard deviations in the RA/Dec residuals of sources deemed 
to have insignificant motion (< 45 milli-arcseconds) between 
epochs. The proper motion of WISEP J075003. 84+272544.8 
is /J-acosS = —732 ± 17mas yr -1 , fig = —194 ± 17mas yr -1 . 
By comparison, Kirkpatrick et al. (2011) used astrometric 
fits to multiple WISE observations to derive a proper motion 
(AWos,5 = -869±424mas yr" 1 , (i s = -1107±438mas yr -1 ) 
with much larger uncertainties. Their value of \i a cos s is con- 
sistent with our new measurement, however their value of 
Us is too large at the level of ~2er. Although there also hap- 
pens to be a nearby Hipparcos star (Hip 38228) , it has a low 
proper motion (n aC osS = — 8mas yr -1 , fig = — lOmas yr -1 ) 
and the T dwarf is not a common proper motion companion 
since its motion differs at a level >28<r. 
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Source Y a J H K Wl(snr) or [3.6] IV2(snr) or [4.5] W3(snr) W4 

WISE 17.75(1. 5) ft 14.76±0.09(11.8) 12.21(0.8) b - 

UKIDSS LAS 19.51+0.14 18.76±0.12 

VISTA VIKING 19.69±0.05 18.71±0.05 

TNG 19.75±0.22 19.14±0.20 

Synthetic Estimate (18.96±0.15) (19.89±0.33) 

Spitzcr 16.77±0.03 14.71±0.01 

a Photometry is on the MKO system except for the TNG Y filter (see text). 
b 95% confidence brightness upper limit. 



Table 3. Photometric magnitudes of WISEP J1423+0116. 



Y — J J - H H - K J - K W1-W2 


W2-W3 


J — W2 


H - W2 


[3.6] - [4.5] 


H - [4.5] 


0.98±0.07 -0.38±0.23 -0.93±0.36 a -1.27+0.34" >2.77 ft 


s£2.55 6 


3.95±0.10 c 


4.38±0.22 


2.06±0.03 


4.43±0.20 


a Synthetic photometry (see text). 












b 95% confidence limit. 












c Using the higher signal-to-noise VISTA J-band magnitude. 













Table 4. Photometric colours of WISEP J1423+0116. 



7 CONFIRMING BINARITY 

To estimate the probability that WISEP J1423+0116 and 
BD+01 2920 may be a line-of-sight association as opposed 
to a genuine binary, we have performed a statistical analysis 
to estimate the expec ted number of chan c e alignm ents in our 
search. We used the iBurningham et al.l l)2010bl ) luminosity 
function constraints to estimate the number of T6-9 dwarfs 
expected in the WISE first data release. In a sample with 
W2 <15 (akin to our W2 signal-to-noise requirement) we ex- 
pect to detect T7±l and T9 dwarfs out to distances (D max ) 
of ~25 and ~15pc respectively, in the 57% sky coverage of 
the WISE first release. We adjusted the Burningham lumi- 
nosity function to add back in the correction made for un- 
resolved binarity (3-45%), since this removed T dwarfs from 
their magnitude limited samples, and estimate an expected 
28-251 T7±l dwarfs and 26-111 T9 dwarfs in the WISE se- 
lection using this luminosity function. We then summed the 
volume in which line-of-sight associated stars may be found 
using a set of cones (one per T dwarf) each with its apex 
at the observer and a T dwarf in the centre of it's base (us- 
ing a base radius of 10,000AU to match our search criteria). 
Since the number of T dwarfs is proportional to D 3 and 
the volume of a cone is proportional to D (where D is dis- 
tance), the average cone volume will be | of the maximum 
cone volume hAD max (where A is the base area of a cone 
7T x 10, 000 AlT). The total cone volume for T6-9 dwarfs was 
thus estimated to be 2.0-14.6pc 3 . 

The luminosity function of lReid et alj (|2007l ) for the 8pc 
and 20pc samples leads to a stellar density of 0.062-0.076 
stars pc -3 , and we thus expect 0.12-1.11 light-of-sight as- 
sociations between stars and T dwarfs in our WISE selec- 
tion. Amongst our five candidates we find that one of them 
(WISEP J075003. 84+272544.8) is in fact a line-of-sight asso- 
ciation with a lack of common proper motion. This is consis- 
tent with our estimates above. An additional candidate was 
identified without proper motion, though the above statis- 
tic does not provide any further indications on the likelihood 
that this candidate may be genuine. 



We must also assess these common proper motion sys- 
tems for the chance that common proper motions are aligned 
by random chance. Using the proper motion and direc- 
tion of WISEP J1423+0116 we estimated this probabil- 
ity using a Hipparcos sample, downloading the proper mo- 
tions of Hipparcos stars within 45 degrees of the WISEP 
J1423+0116/BD+01 2920 pair, and with distances from 
10-40pc (the photometric distance range of a T8±l dwarf 
with J ~18.7 allowing for possible unresolved binarity). We 
counted the fraction of stars with proper motion within 
55mas yr _1 (la) of the T dwarf motion, and thus esti- 
mate a chance of 1.3% that this high proper motion pair 
could be common proper motion by random chance. We 
therefore expect no more than 0.0015-0.014 false positive 
common proper motion systems in the search that we have 
made, and conclude that all three of the common proper 
motion systems that we identified are genuine binaries. This 
includes the two previously reported systems and the associ- 
ation between WISEP J1423+0116 and BD+01 2920, which 
becomes the binary system BD+01 2920AB. 



8 PROPERTIES OF BD+01 2920A 

A search of the literature reveals multiple studies of the pri- 
mary star BD+01 2920A. It is a nearby high proper motion 
Gl dwarf (0.9Mq) with thin disk kinematics. There is no 
evidence of any debris disk or low-mass companions (includ- 
ing giant planets), and it has low activity. BD+01 2920A is 
a mildly metal poor star, with a metallicity in the metal 
poor tail of the disk distribution rather than in the halo 
regime. With one exception pr evious estimates of [F e/H] are 
in the range -0.38+0.06 (only [Eebreton et al.lll999l . gives a 
slightly higher metallicity of -0.20 dex). The range of age 
constraints covers 2. 3-14.4 Gyr. This d iffers slightly from 
the range quoted by lLawler et al. | 1120091) who give a lower 
limit of 3. 5 Gyr. The difference i s due t o the estimate of 2.27 
Gyr from Ido Nascimento et al.l (|2010l ). and we here adopt 
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BD+01 2920A (Hip 70319) 



R.A. (J2000) 


14 23 15.285 


Dec (J2000) 


+01 14 29.65 


P^q cos Dec 


223.8 ± 0.4 mas/yr 


PM Dec 


-477.4 ± 0.4 mas/yr 


Spectral type/class 


G1V 


7T 


58.2±0.5 mas C 1 ) 


Distance 


17.2±0.2 pc 


m - M 


1.18+0.03 


V r 


19.6 ± 0.3 km/s' 2 4 J 


Space motion 


UVW = 22, 15, 39 < 5-8 ' 


Population 


Thin disk ( 9 > 10 ) 


Teff 


5750 ± 100 K (3,7,8,11-21) 


log 9 


4.45 ± 0.05 dex (3,7,12,16-18,20-24) 


Mass 


0.87 ± 0.07 M (3.13,23) 


[Fe/H] 


-0.38 ± 0.06 dex (3,6-8,10,12,14-18,20-29) 


Age 


2 3 - 14 4 Gyr (3, 5, 8, 10,12,13, 15, 23, 26, 30-32) 


vsin i 


1-2 km/s (3.12,13) 


Activity 


Low activity star ( 33 ) 


Debris disk 


None < 15 ) 


Close in companions 


No >70-75M Jup at 20-250AU < 34 > 




No giant planets 




(>100m/s) at < 5AU ( 24 .35-39) 


1 van Leeuwen (2007). 


2 Gontcharov (2006), 3 Valenti & Fischer (2005) 



4 Latham et al. (2002) , 5 Holmbere et al. (2009) , 6 Ramirez et al. (2007) 
7 Mishenina et al. (2004) . 8 Nordstrom et al. (20041 . 9 Borkova fc Marsakov (2004) 
10 Ibukivama fc Arimoto (2002) . 11 Casagrande et al. (2010) . 12 Takeda et al. (2010) 
13 do Nascimento et al. (2010) . 14 Holmbere et al. (2007) , 15 Lawler et al. (2009) 
16 Luck fc Heiter (2006) . 17 Shi et al. (2004) . 18 Mashonkina et al. (2007) 
19 Kovtvukh et al. (2003) . 20 Giridhar fc Goswami (2002) . 21 Cavrel de Strobel et al. (2001) 
22 Wu et al. (2011) . 23 Takeda et al. (2007b) , 24 Heiter fc Luck (2003) 

25 Mashonkina fc Gehren (2001) , 26 Rocha-Pinto fc Maciel (1998) . 27 Karatas et al. (2005) 

28 Borkova fc Marsakov (2005) , 29 Haywood (2001) , 30 Wright et al. (2004) 

31 Barry (1988) . 32 Takeda et al. (2007a) . 33 Hall et al. (2007) 34 Carson et al. (2009) 

35 Halbwachs et al. (2003) . 36 Halbwachs et al. (2000) . 37 Cumming et al. (1999) 

38 Endl et al. (2002) . 39 Nidever et al. (2002) 



Table 5. Properties of BD+01 2920A (Hip 70319). 



an inclusive age range. The properties of BD+01 2920A are 
summarised in Tableland references therein. 

9 PROPERTIES OF BD+01 2920B 
9.1 Bolometric flux 

We estimate the bolometric flux (Fbol) of the new T8p com- 
panion BD +01 2920B followin g a sim ilar method to that 
outlined in iBurningham et all |2009a), by combining our 
YJHK spectrum (flux calibrated by our follow-up photom- 
etry) with model spectra (to allow us to estimate the flux 
contributions from regions outside our near-infrared spec- 
tral coverage). We have scaled the A < 1.05(im region of 
the models to match the flux level in our YJHK spectrum, 
and we have used the Spitzer 3.6um and 4.5(J.m photome- 
try to scale the A = 2.5 — 3.95[im and A > 3.95um regions 
respectively (the transmission profiles of the Spitzer filters 
cross at 3.95 urn at a transmission level of ~1%). To avoid 
biasing our derived flux estimate with our choice of model, 
we have produced estimates using Solar and [M/H]= —0.3 
metallicity BT Settl models (|Allard et al.ll201ll ) that bracket 
the likely range of gravities and T e ff for our target (logp = 



4.50 - 5.5; T eS = 500 - 900K). We take the scatter in the es- 
timates resulting from different model choices as a reflection 
of the systematic uncertainty introduced by the atmospheric 
models. We have used a Monte Carlo method to determine 
the uncertainty in each estimate due to the noise in the 
photometry used for scaling the models and the noise in our 
GNIRS spectrum. Our final estimate of Fbol is the median 
of our estimates using different model extensions, whilst the 
uncertainty is the sum in quadrature of the systematic un- 
certainty and the mean random uncertainty. This results in 
a determination of F bol = 1.61 ± 0.17 x 10~ 16 Wm~ 2 . 

9.2 Luminosity, mass, radius, and effective 
temperature 

The luminosity of BD+01 2920B was derived from the bolo- 
metric flux and the distance. The on-sky separation of the 
BD+01 2920AB components leads to a tangential separation 
of 0.01 pc. This is negligible compared with the uncertainty 
in the parallax distance of the primary (±0.2 pc) and we can 
therefore assume that the T dwarf is at the same distance as 
BD+01 2920A (17.2+0.2 pc). Taking the uncertainties as- 
sociated with the bolometric flux and distance into account 
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BD+01 2920B (WISEP 1423+0116) 



R.A. (J2000) 
Dec (J2000) 

P^a cos Dec 

PM Dec 
Spectral type 
Separation 

Fbol 

m — M 

My 

Mj 

M H 

M K 

M 3 .e 

M4.5 

log L/Lq 

[Fe/H] 

Mass 

Radius 
logs 



14 23 20.86 

+01 16 38.1 

261 ± 56 mas yr _1 

-444 ± 52 mas yr~ 

T8p 

153 arcsecs 
2630 AU a 
(1.61±0.17)xl0- 16 
1.18+0.02° 
18.51+0.04 a 
17.53+0.05 a 
17.96+0.20 a 
18.71+0.33 a - c 
15.59+0.04 a 
13.53+0.02 a 
-5.83 ± 0.05 a 
-0.38 ± 0.06 dex d 
0.019-0.047 M e 
20-50 M Jup e 
0.080-0.099 R e 
0.80-0.99 R Jup e 
4.68-5.30 dex e 
680+55 K / 



Wm 



a Inferring a distance of 17.2±0.2pc from BD+01 2920A. 
b Integrating measured flux from 1.0 — 2.4um and adding a 
theoretical correction at other wavelengths (see text). 
c Synthetic photometry used (see text). 
d Inferred from BD+01 2920A. 

e Constraints derived from structure models as a function of 
luminosity for ages 210 Gyr. 

* Derived from the luminosity and radius constraints. 
Table 6. Properties of BD+01 2920B (WISEP 1423+0116). 



leads to the determination of L bot = 5.69 ± 0.60 x 10" 20 W, 
or log L/Z/Q = -5.83 ± 0.05. To determine the T e g of 
BD+01 2920B we estimate its radius using the COND evo- 
lutionary models jBaraffe et al. 2003). These models repro- 
duce the main trends of observed methane dwarfs in near- 
infrared color-magnitude diagrams, though are only avail- 
able for solar metallicity. 

We used linear interpolation between the model 
isochrones to estimate a range of possible mass, radii and 
surface gravity for BD+01 2920B consistent with an age 
range of ~2-10 Gyr. Accounting for the uncertainties in the 
measured luminosity we obtained a mass range of 0.019- 
0.047 Mq (20 - 50 M Jup ), a radius range 0.080-0.099 R 
(0.80-0.99 Kjup), and a surface gravity range of loggr = 
4.68 — 5.30. The corresponding temperature (from luminos- 
ity and radius) is T e g = 680 + 55K. A summary of the prop- 
erties of BD+01 2920B is given in Tabled 

Observations of transiting very low-mass stars an d 

~ 2005bllal; 



brown dwarf s with mass>20 Mj v „ llPont et al 



Deleuil et all |200Sl: iBouchv et al.l 1201 lbl : I Anderson et all 



201ll ; |johnson et al.ll201ll ) are all "cons istent with the COND 



mass-radius model data (see fig 10 of lBouchv et al.ll2011al ). 
These systems are solar metallicity to within the uncertain- 
ties, though some of these uncertainties are significant. The 
effects of meta ll icity on sub-stellar radii are a little unclear. 
I Burrows et al.l (|201ll ) present evolutionary models with a 
spread in radius (at a given mass and age) of ~10-25%, with 



higher-metallicity (higher-cloud-thickness) atmospheres giv- 
ing larger radii. However, a comparison between KOI-423b 
and CoRoT-3b suggests the converse trend. KOI-423b orbits 
a metal poor star ([Fe/H] =— 0.29+0.10) and is a relatively 
large (1.22t.o'\oRjup) brown dwarf (18 Mj up ), whereas 
CoRoT-3 is a solar metallicity star ([Fe/H]=-0.02±0.06) 
hosting a smaller (1.01±0.07Rj up ) brown dwarf (22Mj up ). 
Given this uncertainty in the radius-metallicity trend, we do 
not attempt to make a metallicity correction to our radius 
estimate. We note however that our COND radius constraint 
already includes an uncertainty at the level of 25%, compa- 
rable with the s ize of the theoretical trends suggested by the 
I Burrows et al.l (|201ll ) models for a metallicity difference of 
1.0 dex. 

As an additional caveat we note that our T e g determi- 
nation relies on an assumption that the object is single, and 
not an unresolved binary. Unresolved binarity would lead to 
lower T e ff for each unresolved component. If BD+01 2920B 
is an equal luminosity unresolved binary the T e g of each 

component would be a factor ~0.8 less (^ 3 with similar 
radii for the components). For unequal luminosity com- 
ponents the brighter component T e g would be closer to 
680K with the fainter one <540K. Observations suggest (e.g. 
iBurgasser et~ai1 120051 ) that the binary fraction of brown 
dwarfs (resolved at ~0.1 arcsec resolution) in widely sep- 
arated stellar - brown dwarf multiples is notably higher 
(45+15 per cent) than that of field brown dwarfs (18+7 
per cent), and unresolved binarie s can also have separation 
closer than 0.1 arcseconds (see iBurningham et all l2009al . 
and references therein). In Figure[3]we show BD+01 2920B 
in absolute magnitude (Mj,h,k) spectral type diagrams, 
along with the known population of late T and Y dwarfs 
(see caption). The iC-band suppression is evident in the Mk 
plot, though we also note that there is no clear indication 
that the object is an unresolved binary (e.g. with compo- 
nents of near equal brightness) in the Mj and Mr plots. 
We cannot with high confidence however, rule out the pos- 
sibility that BD+01 2920B is an unresolved binary. 



10 A METAL POOR BENCHMARK T8 
DWARF 

We now assess some implications of this benchmark system 
under the assumption that it is a single object, and through 
comparison of its spectrum and colours to theoretical pre- 
dictions and other ultra-cool objects. In Figure [4] we com- 
pare our flux calibrated GNIRS spectrum of BD+01 2920B 
and warm- Spit zer photometry t o mildly metal poor BT 
Settl models (jAllard et al.l [2010) for our derived proper- 
ties, each scaled to their corresponding radii and the known 
distance to the primary star. The BT Settl atmospheric 
model grid spans the cool stellar to su bstellar temperature 
regime using the BT2 water line-list (Barber et al. 2006) 
and the reference Solar abundances o f lAsplund et al.l l|2009l) . 
In these models dust formation, cloud behaviour and verti- 
cal mixing are parameterised with re ference to the 2D ra - 
diation hydrodynamic simulations of iFrevtag et al.l l|2010t ). 
The K— band suppression that is predicted by the models 
for high-gravity and metal poorer brown dwarfs is seen in 
our GNIRS spectrum of this benchmark T dwarf, although 
the models predict this effect to be stronger than is seen in 
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• BD+01 2920B 

• Single objects (32) 
Known binories (2) 
Possible binaries (1) 




IR Spectral type 

Figure 3. Absolute magnitude (Mj } h,k) spectral type plots 
showing BD+01 2920B (red symbol) amongst the existing popu- 
lation of mid-late T d warfs and Y dw arfs. All the T dwarfs with 
spectral types ^T8.5 jBurningham e t al. 2008|) have distances 
from parallax. The T9.5 and Y dwarfs all have spectroscopic dis- 
tances (except for WIS EP J1541-2250 which has a parallax in 
iKirkpatrick et al.|[2uTiT) , estimated via T c ff and log g constraints 
from model fits to near-infrared spectroscopy. 



this case. Similarly poor matches to the observed K— band 
spec troscopy have been seen in othe r benchmark systems 
(e.g. IBurningham et al.ll2009bl . l2011bl ) , so it is reasonable to 
interpret this as a deficiency in the model atmospheres, al- 
though its origin is not yet understood. It is noteworthy that 
the Y— band spectrum which has also been proposed as di- 
agnostic of metallicity variations (IBurgasser et alj|2006ar ) is 
also poorly matched by the models. The model atmospheres 
provide good matches to the flux in the J— and H— bands, 
despite known deficiencies in the methane line lists in these 
regions. The model predictions for the [3.6] and [4.5] fluxes 
are consistent with those that we have observed, and it inter- 
esting to note that in the [3.6] band, the 620K, loggr = 4.7, 
model fits the data best, whilst in the [4.5] band the two 
warmer models provide the best match. 

In Figure [5] we have plotted synthetic colours in 



RD+01 ?9?0B 

BT Settl 620K, log g = ' 

BT Settl 580K, log g = t 

BT Settl 740K, log g = 5.3 




Figure 4. The YJHK spectrum of BD+01 2920B and mean 
fluxes inferred from the Spitzer photometry compared to model 
spectra that straddle the properties estimated in Section 19.21 The 
straight coloured lines indicate the mean fluxes of model spectra 
in the IRAC channel 1 and 2 photometric bands and are plotted 
to allow comparison with the mean flux from the target (straight 
black lines) through the same filters. 



H - K and H - [4.5] for the B T Settl mode l s along 
with colours of late T dwarfs from iLeggett et al.l l|2010b ) 
with MKO and IRAC photometry, and BD+01 2920B. 
It can be seen that BD+01 2920B lies in a similar re- 
gion as other suspected mildly metal-poor T7.5/8 dwarfs 
2M ASS0939 and SDSS1416B which have Teff 5 00-700K 
(e.g. lBurningham et al . 2010a: Leggett et al 1l2010ah . The ef- 
fect of the poor match between the models and the data in 
the if— band is highlighted by the non-coincidence of the 
models and the observations in this colour space for the T 
dwarfs with the reddest H — [4.5] colours. However, the mod- 
els correctly predict the colours for the young benchmark 
Ross 458C. To provide an alternative comparison between 
the models and the data we have shifted the model colours 
such that they match the observed colours for BD+01 2920B 
for the parameters derived in Section [9] Figure [6] com- 
pares these adjusted model colours to the same T dwarfs 
sh own in Figure [ 5] This plot broadly reproduces the result 
of ILeggett et al.l pHObh . where it was noted that the ma- 
jority of the coolest T dwarfs appear to have low-gravity 
and/or high- metallicity, suggesting that the sample is domi- 
nated by young low- mass brown dwarfs (ages ~ lGyr). How- 
ever, we note that the sources with bluer H — [4.5] lie well 
below the adjusted model tracks in Figure [6l including the 
young benchmark Ross 4 58C (for which log;; = 4.0 — 4.7; 
IBurningham et all l2011al ) , which highlights that a simple 
offset correction to the models is not sufficient to allow 
the properties of the T dwarf population to be reliably as- 
sessed through reference to these model colours. The similar 
temperature of these two benchmarks (T c ff ~ 700K), but 
wide separation in Figures and [5] highlights the important 
roles that gravity and metallicity play in determining the 
H — K/H — [4.5] colours for cool T dwarfs. 
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H - [4.5] 



Figure 5. Near- to mid-infrared colours of cool T dwarfs com- 
pared to those of the BT Settl model colours. BD+01 2920B is 
indicated with a cyan star symbol, whilst the young benchmark 
Ross 458C is indicated in red. 




2.5 3.0 3.5 4.0 4.5 5.0 

H - [4.5] 



Figure 6. Near- to mid-infrared colours of cool T dwarfs com- 
pared to those of the BT Settl model colours anchored to our 
estimated properties for BD+01 2920B. BD+01 2920B is indi- 
cated with a cyan star symbol, whilst the young benchmark Ross 
458C is indicated in red. 



11 CONCLUSIONS AND FUTURE WORK 

Our cross-match between the WISE first data release and 
the Hipparcos and Gliese catalogues has resulted in the dis- 
covery of a new late T binary companion (BD+012920B) 
and the re-discovery of two previously known systems. WISE 
(in combination with UKIDSS and VISTA) is thus effec- 
tively probing an increased volume of very low temperature 
parameter-space for benchmark companions. There are also 
significant advantages that the primary star BD+012920A 
is a nearby G dwarf rather than one of the more numerous 
M dwarfs in the solar neighbourhood, and its lower metal- 
licity provides a crucial test for the effects of reduced metal 
content on models atmospheres. The metallicities and abun- 
dances of bright Sun-like stars can be studied with much 
more confidence and in much greater detail than those of M 
dwarfs, and late T dwarfs in such systems offer the opportu- 
nity not only to test cool brown dwarf atmosphere physics, 
but also to potentially study brown dwarf abundances. 

In the near future, high resolution imaging (e.g. adap- 



tive optics) observations of BD+012920B will be important 
to constrain multiplicity on a ~0.1 arcsec (~1.7AU) sepa- 
ration scale. A close binary would be able to provide future 
dynamical masses (since the orbital period would be just 
a few years). Higher resolution spectroscopy would also be 
able to assess multiplicity at closer separation, and confir- 
mation of a single object nature would validate the approach 
taken here to determination the physical properties of this 
benchmark object. The existing constraints on the physical 
properties of BD+012920B will be improved as we develop 
a better understanding of how brown dwarf radii depend 
on composition. This will be aided by an increasing number 
of transiting brown dwarfs from Kepler and other t ransit 
surveys (e.g. iBorucki et alll201ll ; IPinfield et aT1l2005h . and 
improved metallicity measurements for this sample. 

We can expect additional late T benchmarks in the fu- 
ture all-sky WISE data release, and a more encompassing 
search of WISE, UKIDSS and VISTA (including at wider 
angular separations) should yield an expanded population 
of benchmarks across the full T dwarf T e a range. As greater 
survey volumes are searched for benchmark brown dwarfs 
we can also expect to identify systems with more accu- 
rately known ages. Evolved subgiants for example are less 
numerous that their main sequence counterparts, but evolu- 
tionary model compariso ns can provide mor e accurate age 
constraints (e.g. ±10%; [Pmfield et al.l l200fj ). And as the 
range of well measured benchmarks expands into greater 
parameter-space we will have the opportunity to compre- 
hensively test the atmosphere models by directly mapping 
the benchmark population's spectral variations/trends onto 
a grid of tightly constrained physical properties. 
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